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Abstract 

We investigate the charmed baryon mass spectrum using the relativistic heavy quark action on 
2+1 flavor PACS-CS configurations previously generated on 32 3 x 64 lattice. The dynamical up- 
down and strange quark masses are tuned to their physical values, reweighted from those employed 
in the configuration generation. At the physical point, the inverse lattice spacing determined from 
the Q baryon mass gives a" 1 = 2.194(10) GeV, and thus the spatial extent becomes L = 32a = 
2.88(1) fm. Our results for the charmed baryon masses are consistent with experimental values, 
except for the mass of H cc , which has been measured by only one experimental group so far and 
has not been confirmed yet by others. In addition, we report values of other doubly and triply 
charmed baryon masses, which have never been measured experimentally. 
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I. INTRODUCTION 



Recently, a lot of new experimental results are reported on charmed baryons [1]. BaBar 
and Belle give very accurate results based on their precise analysis. In addition, new ex- 
periments such as J-PARC, PANDA, LHCb, and Belle II are coming and expected to give 
further informations for charmed baryons. 

Mass spectrum of singly charmed baryons has been determined experimentally in high 
accuracy. Experimental status for masses of the ground state is evaluated as three or four- 
star by the particle data group. The excited states are also investigated fairly well. 

In contrast to singly charmed baryons, experimental data for doubly and triply charmed 
baryons are not well established. A candidate for the doubly charmed baryon, S cc , has been 
reported only by the SELEX Collaboration [2J, while S cc has not been confirmed yet by other 
experimental groups such as BaBar [3| and Belle Q. Further experimental and theoretical 
confirmations are required to establish S cc . Other doubly charmed baryons and any triply 
charmed baryons have not been observed by experiments. In this situation, theoretical 
predictions for the doubly and triply charmed baryon masses may give useful information 
for the experimental discovery of these states. 

Charmed baryon spectrum has been mainly investigated using gauge configurations gen- 
erated with 2+1 flavors dynamical staggered quarks [5|-l8|. In this case, a choice for light 
valence quarks requires a special care: One may take other fermion formulations for the va- 
lence light quarks (the mixed action) to avoid a problem due to a complicated flavor structure 
of the staggered quarks. Alternatively one may construct charmed baryon operators with 
the valence naive quark and then rewrite the correlation function of these operators in terms 
of the staggered propagators [9|. Both approaches violate the unitarity of the theory at finite 
lattice spacing, in addition to the rooting problem of dynamical staggered quarks. There- 
fore it is necessary to check their results using other combinations of dynamical and valence 
quarks which maintain the unitarity at non-zero lattice spacing. Another issue in the ex- 
isting calculations on dynamical staggered configurations is that their chiral extrapolations 
suffer from large higher order corrections. NLO SU(2) heavy baryon chiral perturbation 
theory is employed to extrapolate their data to the physical point from 220 - 290 MeV 
pion masses, but the result shows a bad convergence even at m v = 220 MeV. A calculation 
directly at the physical quark masses without chiral extrapolation is the best way to remove 
this uncertainty. 

There are a few investigations with other fermion formulations. S. Diirr et al have calcu- 
lated charmed omega baryon masses using smeared improved Wilson and Brillouin fermions 
for valence strange and charm quarks |10f] . Though their calculation is performed only at 
one quark mass (m n = 280 MeV) and one lattice spacing (a = 0.07 fm) on the Nf = 2 
0(a)-improved Wilson quark ensemble, they obtained a result which is consistent with the 
experimental value for the singly charmed omega baryon (Q c ) mass. 

The ETMC has studied charmed baryons on ensembles generated with Nf = 2 dynamical 
twisted-mass quarks, employing the same twisted-mass fermion for degenerate up and down 
valence quarks and Osterwalder-Seiler fermions for strange and charm valence quarks (llj . 
For a doubly charmed baryon, they have found m= cc = 3.513(23) (14) GeV. This is the only 
result that is consistent with the SELEX experimental value m|^ c LEX = 3.519(1) GeV, while 
other lattice QCD results deviate from this value. Reasons for this disagreement among 
lattice QCD results must be understood and should be eventually resolved. One possible 
source for systematic uncertainties in the ETMC calculation is a lattice artifact caused by 
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the heavy charm quark mass at their lattice spacings, a = 0.09 — 0.06 fm. Indeed their results 
for charmed baryon masses, especially for rriQ ccc , do not show clear scaling behaviors. To 
reduce this uncertainty, one must employ a heavy quark action that handles mass dependent 
lattice artifacts in the formulation, such as the Fermilab action [12], the relativistic heavy 
quark action 13J, |l4| , or highly improved actions. Chiral extrapolation of ETMC data from 
= 260 MeV using the NLO heavy baryon chiral perturbation theory is another source 
of systematic uncertainties, as in the case of staggered quarks. 



In Ref. [15], the present authors have shown that the charm quark mass corrections are 
under control at a -1 = 2.194(10) GeV by adopting the relativistic heavy quark action of 
Ref. [13]. It removes the leading cutoff errors of 0((mQa) n ) and the next to leading effects of 
0((rriQa) n (aAQCD)) for arbitrary order n by tuning finite number of parameters. Employing 
this action for the charm quark, we have investigated properties of mesons involving charm 
quarks with the 2+1 dynamical flavor PACS-CS configurations on 32 3 x 64 lattice [l6[ 
reweighted to the physical point for up, down and strange quark masses. We have found 
our results for charmed meson masses are consistent with experiment at a percent level, and 
so are those for the decay constants with a few percent accuracy, though our results are 
obtained at a single lattice spacing. 

Encouraged by this result, we have extended our investigation to the charmed baryon 
sector and report the result in this paper. One of the big advantage in our investigation 
over previous calculations is that the chiral extrapolation is no more necessary, since our 
calculations are performed at the physical point. We are free from the convergence problem 
of the heavy baryon chiral perturbation theory. We first compare our masses of singly 
charmed baryons with the corresponding experimental values, to check if our method works 
also for the baryon sector. We then evaluate the doubly and triply charmed baryon spectra 
as our predictions. A part of this work has been reported in Ref. (l7| . 

This paper is organized as follows. Section II explains our method and simulation pa- 
rameters. Section III describes our results for the singly charmed baryon spectrum, and 
comparison with experiments. In Sec. IV, we present our results for doubly and triply 
charmed baryon masses. A conclusion is given in Sec. V. 



II. SET UP 

Our investigation is based on a set of 2 + 1 flavor dynamical lattice QCD configurations 



generated by the PACS-CS Collaboration 16J on a 32 3 x 64 lattice using the nonpertur- 
batively_0(a)-improved Wilson quark action with Cg^ = 1.715 [l8[ and the Iwasaki gauge 



action [19J at (3 = 1.90. The aggregate of 2000 MD time units were generated at the hopping 
parameter given by k°) = (0.13778500,0.13660000), and 80 configurations separated 
by 25 MD time units were selected for our calculations. We then reweight those configura- 
tions to the physical point given by (n u d, K s) = (0.13779625,0.13663375). The reweighting 
shifts the masses of tt and K mesons from = 152(6) MeV and m# = 509(2) MeV to 
= 135(6) MeV and m# = 498(2) MeV, with the cutoff at the physical point estimated 
to be cr 1 = 2.194(10) GeV from the Vt baryon mass. Our parameters and statistics at the 
physical point are given in Table [H 

The relativistic heavy quark action [l3[ is designed to reduce cutoff errors of 0((m<ga) n ) 
with arbitrary order n to 0(f (rriQa)(aAQCD) 2 ) , once all of the parameters in the action are 
determined nonperturbatively, where /(mga) is an analytic function around the massless 
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# conf MD time 



1.90 0.13779625 0.13663375 80 2000 

TABLE I: Simulation parameters. MD time is defined as the number of trajectories multiplied by 
the trajectory length. 



^charm V_ Vg CB C E 

0.10959947 1.1450511 1.1881607 1.9849139 1.7819512 
TABLE II: Parameters for the relativistic heavy quark action. 



point rrtQa = 0. The action is given by 
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The parameters t s ,cb,ce and v have been tuned in Ref. [15|. It should be noticed that, 
while the perturbative estimate are used for t s ,cb,ce [2(j, the parameter v is determined 
non-perturbatively to reproduce the relativistic dispersion relation for the spin-averaged 15* 
state of the charmonium. The heavy quark hopping parameter kq is set to reproduce the 
experimental value of the mass for the spin-averaged IS state. Our parameters for the 
relativistic heavy quark action are summarized in Table [TTJ 

We employ the relativistic operators for charmed baryon, simply because the relativistic 
heavy quark action is employed in our calculation. Charmed baryons can be classified under 
4 x 4 x 4 = 20 + 20i + 20 2 + 4. In addition to J = 3/2 decuplet-type 20-plet, there are 
J = 1/2 octet-type 20-plet and 4-plet. 

J = 1/2 octet-type baryon operators are given by 



O f a 9h ( 



x 



e^((q a f (x)y C l5 q b g (x))q c ha (x) 
C = 7472, 



(II.3) 
(II.4) 



where /, g, h are quark flavors and a, b, c are quark colors, a = 1, 2 labels the ^-component 
of the spin. The S-type and A-type are distinguished as 



S-type : 
A— type : 



_|_ QWf 

Q[fh]g _ Q[gh]f _ 20 

7$ ' 



(II.5) 
(II.6) 



where 0^ f9]h = O fgh - 9fh . 
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TABLE III: Our results for singly charmed baryon masses. Experimental values from PDG [l( are 
also listed. 



The decuplet-type J = 3/2 baryon operators are expressed as, 
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r± = (71 =R72)/2,r = 73. 


(11.11) 



Two-point functions for these charmed baryon operators are calculated with exponentially 
smeared sources and a local sink. The smearing function is given by \l/(r) = Aexp(-Br) at 
r ^ and (0) = 1. We set A = 1.2, B = 0.07 for the ud quark, A = 1.2, B = 0.18 for the 
strange quark, and A = 1.2, B = 0.55 for the charm quark. The number of source points is 8 
per configuration and polarization states are averaged over to reduce statistical fluctuations. 
Statistical errors are analyzed by the jackknife method with a bin size of 100 MD time units 
(4 configurations), as in the light quark sector (l6T |. We extract charmed baryon masses 
by fitting two-point functions with single exponential forms. Figure dHS] show our effective 
masses. We take the fitting interval for charmed baryons as far as possible from the origin, 
[tmin,t m ax] = [10, 15], to avoid possible contaminations from excited states. Our results are 
compiled in Table II I II and IIVI 



III. SINGLY CHARMED BARYON SPECTRUM 



Our results for the singly charmed baryon spectrum at the physical point are summarized 
in Fig. HI All our values for the charmed baryon masses are predictions from lattice QCD, 
since the physical charm quark mass has already been fixed by the mass for the spin-averaged 
IS charmonium state and no other experimental inputs for charmed baryon masses are 
required. 
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FIG. 1: Effective masses of J = 1/2 singly charmed baryons. 



As can be seen from the figure, our prediction for the singly charmed baryon spectrum in 
lattice QCD at a single lattice spacing is in reasonable agreement with the experimental one. 
In Fig. |5l we also compare our value for A c with other results obtained in recent lattice QCD 
simulations using the dynamical staggered quarks [B|-0], and the twisted mass quarks (llj . 
All results are consistent with each other, though the statistical error is larger for our result 
due to the conservative choice of our fitting interval. 

Figure [6] displays mass differences. Our results are consistent with experimental values 
within 2 o uncertainty. These agreements indicate that the decomposition of J — 1/2 E- 
type and A-type baryons, as well as that of J = 3/2 and J = 1/2 charmed baryons, have 
been made successfully in our calculation. 

We note that several systematic errors have not been fully evaluated yet for these results. 
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FIG. 2: Effective masses of J = 3/2 singly charmed baryons. 
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TABLE IV: Our results for doubly and triply charmed baryon masses. An experimental value 
from SELEX [2] is also listed. 



First, finite size effects are not taken into account. Although the NLO heavy baryon chiral 
perturbation theory predicts that finite size effects for charmed baryons are less than 1 
%, higher order terms may give significant contributions. A direct confirmation in lattice 
QCD by comparing spectra among different lattice volumes is desirable. Second, strong 
decays such as S c — > A c ir are not taken into account in our analysis, since S c — > A c tc is 
kinematically prohibited on our lattice volume. Last but not least, our results are obtained 
at a single lattice spacing without continuum extrapolation. Although a naive order counting 
gives a percent level of cutoff effects from 0(a 2 s f{mQa)(aA.QCD), f( m Qa)(aAQCD) 2 ) terms in 
the relativistic heavy quark action, the continuum extrapolation is necessary to remove 
this uncertainty. Additional calculations should be performed in the future to remove all 
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FIG. 3: Effective masses of doubly and triply charmed baryons. 



systematic errors mentioned above. 



IV. DOUBLY AND TRIPLY CHARMED BARYON SPECTRUM 



For doubly and triply charmed baryons, an experimental value has been reported only 
for S cc , though the experimental status is controversial. In the other channels, lattice QCD 
result gives predictions before experimental mass measurements. 

Figure [7] shows our results for the doubly charmed baryons. Our estimate for m= cl . 
clearly deviates from the experimental value by SELEX Collaboration (3] . The difference is 
4cr, as shown in the right figure. Our result for m= cc is consistent with results from other 
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FIG. 4: Our results for the singly charmed baryon spectrum (left panel), and those normalized 
by the experimental values (right panel). 
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lattice QCD calculations except ETMC, as shown in Fig. [SJ This discrepancy is need to be 
understood and should be resolved. 

Similarly, Fig. [9] displays lattice QCD results for the triply charmed baryon from several 
groups. Our prediction agrees with that by others except ETMC. A marginal discrepancy 
is observed between our value and that of Ref. [§] for mcj ccc — 3/2 rrij/^ mass difference (the 
right figure). 

For a more detailed comparison, precise evaluations of all systematic errors are required. 
In particular, the largest source of systematic errors for our calculations is the lattice artifact, 
which should be removed by the continuum extrapolation. 



V. CONCLUSIONS 



We have studied charmed baryon masses in Nf = 2 + 1 dynamical lattice QCD at a 
lattice spacing of a -1 = 2.194(10) GeV. The reweighting technique allows us to perform a 
measurement directly at the physical point. It removes a systematic error associated with 
chiral extrapolations of charmed baryon masses, which had prevented previous lattice QCD 
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FIG. 7: Our results for the doubly charmed baryon spectrum (left panel), and those normalized 
by the experimental value (right panel). 



calculations from predicting precise values for charmed baryon masses. 

Our results for the mass spectrum of singly charmed baryons are consistent with experi- 
ments within 2 a uncertainty. This confirms that we are able to control charm quark mass 
corrections successfully, not only in the meson sector [TEJ but also in the baryon sector, by 
use of the relativistic heavy quark action of Ref . [l3[ . 

We then extract predictions for doubly charmed baryons. Our result for m= cc is consistent 
with values of other lattice QCD calculations employing the dynamical staggered quarks, 
but disagree with the estimation by the ETMC. Moreover, our S cc mass is different from the 
SELEX experimental value, approximately by 85 MeV, which corresponds to 4a. A similar 
deviation between ETMC and us is also observed in the triply charmed baryon mass, mri ccc . 
Precise estimations of all systematic errors, especially the lattice artifacts, are required to 
resolve these discrepancies. 
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